On the ferromagnetism equations with large variations solutions 

Olivier Gues and Franck Sueur^ 
Abstract 

We exhibit some large variations solutions of the Landau-Lifschitz equations as the ex- 
' change coefficient tends to zero. These solutions are described by some asymptotic ex- 

. pansions which involve some internals layers by means of some large amplitude fluctuations 

in a neighborhood of width ^ e of an hypersurface contained in the domain. Despite the 
nonlinear behaviour of these layers we manage to justify locally in time these asymptotic 
■ expansions. 
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O '. 1 Introduction 
cn 

Ferromagnetic materials can attain a large magnetization under the action of a small applied 
PLh ' magnetic field. To explain this phenomenon, in 1907, Weiss suggested that a spontaneous mag- 

■ netization occurs. In 1928 Heisenberg explained the spontaneous magnetization postulated by 

fH , Weiss in terms of the exchange energy. In 1935 Landau and Lifschitz (cf. proposed a quan- 

I titative theory, now known as micromagnetics. For a piece of ferromagnet-which is supposed to 

be a regular bounded open set Q. in R'^ with a smooth boundary, and locally on one side of F- 
the magnetic state at a point x S ^2 at time t is given by a vector n(t, x) G which belongs to 
the unit sphere of M^, called the magnetic moment. The Landau-Lifschitz equations read: 

d^u^ = -u" A {rL{u^) + e^Au") - A {u^ A {rL{u^) + e^An^)) in (1.1) 

(N 

\ where e > is the exchange coefficient. We denote 'H(u) := H\q^ G L^(il;M'^) where the magnetic 

' field H G L^(M^;R^), is the unique solution of the following elliptic problem, 

\D ' 
O 



curl = in , (1.2) 
div (F + u) = in M\ 



^ , where u means the extension of u by outside of the set The equations (|1.1|) are supple- 

^ ' mented by the homogeneous Neumann boundary condition: 



dnu" = in F, (1.3) 
where n is the unitary outward normal at the boundary F, and by an initial condition: 

4=0 = ^0- (1-4) 
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The solution must also satisfy the constraint 



\u%t,x)\ = l, ioTxen,t>0 (1.5) 

which is obviously propagated from the initial data as soon as it is satisfied at t = 0. 

In this paper we study the asymptotic behaviour of the solutions of the Landau-Lifschitz 
equations (|l.lj) - (|1.3j) - (|1.4jl as the exchange coefficient e tends to zero. From a formal point of 
view, when e = 0, the system (|l.l|) - p.;-{j) - p.4|) becomes 

' dtu^ = n° A n{u^) - u° A (n° A 7^(n°)) in Q 

(1.6) 

where no boundary condition is needed. In the paper |^ it is proved that, for smooth enough 
solutions the system (|1.6|) is a "good approximation" of the full system (|l.lj) - (|1.3j) - (|1.4|) in the 
sense that the solution of ()1.6() is indeed limit in L^([0, T] x Q) of solutions of (jl.lj) - 
(|1.3|1 - (|1.4() as e ^ 0. However, this result holds under the assumption that belongs to the 
space C[[0,T], H^{Q)^ where H^{Q) is the usual Sobolev space. In particular this assumption 
excludes the case where is discontinuous across an hypersurface contained in Q and it was 
one motivation behind this paper to treat that case. 

First of all, let us observe that the system (|1.6|) actually admits discontinuous solutions. To 
simplify, we will restrict the analysis to piecewise smooth solutions. We assume that S is a 
smooth compact hypersurface contained in Q. For < s < cxd call H^{Q) the set of functions 
u G L^(r2) such that u^q^ € H''{Q±) where H^{0.±) is the usual Sobolev space on L^. We endow 
H^{Q) with the norm 

\\u\\hi, ■= \\u\n^\\H'{n-) + 

This definition extends to the case when s = oo: the space H^{Q) is the natural Frechet space. 
We get the following result of global existence of solution of ()1.6() discontinuous through the 
hypersurface S. 

Theorem 1.1. Let s g]|,oo] and uq e H§.{fl). Then there exists a unique vP G (M, if|;(r2)) 
solution of the Cauchy problem (| j.6|) . 

Proof. This result can be easily obtained by following the proof of Proposition 4.1 of I^, with only 
a few adaptations. By the way in the closer setting of semilinear symmetric hyperbolic system, 
it is well known since the works of Metivier ^| that there exist some local piecewise regular 
solutions discontinuous across a smooth hypersurface which is a characteristic hypersurface of 
constant multiplicity for this hyperbolic system. In the present setting, the proof is in fact 
simpler since the hypersurface S is totally characteristic. Moreover thanks to (|1.5() and since 
the operator H satisfies the transmission property, our setting allows to conclude to a global 
existence. □ 

Let us now claim a first theorem about the asymptotic behaviour of the solutions of the 
Landau-Lifschitz equations Hl.l|) - ()1.3() - H1.4|) as the exchange coefficient e tends to zero. 

Theorem 1.2. Let u° G C°°(M, H^{n)) be a solution of (fTTfij) . There exist T > and a family 
of solutions £ C°°([0,T] x 0,), e g]0, 1], of the equation on [0,T] x fi, of the equation 

(| j.,^ on [0,T] X r, such that there exist C > and eq such that for all e G]0,eo], 

- u°\\Lmo,T]xQ) ^ Ce'^. 
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To begin with some comments about Theorem II .21 let us stress that we do not prescribe the 
initial data (|1.4() for the . Thus the traces of the a.t t = are not equal in general to the 
trace of at t = 0. So Theorem ll.2l claims the existence of local in time solutions £ C°°, of 
the equation on Q, of the equation H1.3|) on T, which converge to as e tends to zero in 
L^, with a rate of convergence in £2. 

Indeed, in this paper, we will claim a more accurate result in Theorem I2.1l bv showing that 
the u'^ can be described with a WKB expansion which involves some boundary layers profiles. On 
one hand, a boundary layer appears near the boundary to compensate the lost of the Neumann 
condition from the complete model (|1.1() - (|1.3() - H1.4|) to the limite model ()1.6() (e = 0). Such a 
boundary layer was already studied in paper gj. The amplitude of this boundary is weak and 
its behaviour is linear. On another hand, there are some boundary layers on each side of the 
hypersurface S. Their task is to compensate the lost of transmission conditions across S from 
the complete model (fTTT|) - (fO|) - (|rH) to the limit model (fO)]) {e = 0). 

Remark 1.1. Such an analysis is inspired by the paper jl3j where we show that discontinuous 
solutions of multidimensional semilinear symmetric hyperbolic systems, which are regular outside 
of a smooth hypersurface characteristic of constant multiplicity, are limits, when e — > 0, of 
solutions (M^)jg]o,i] of the system perturbated by a viscosity of size e. In this paper, we adapt 
the method to the ferromagnetism quasi-static model, where in particular the non local operator 
Tl occurs. We point out that for the limit model (e = 0), the hypersurface is totally characteristic. 
As a consequence, the analysis involves only characteristic boundary layers. On the opposite, 
jl3j stresses the occurrence of characteristic and non characteristic boundary layers. It could be 
also possible -as in ^3]- to study the case where the singularity is weaker than a jump of the 
function as a jump of a derivative of the function u^. Then we can take T as big as we want 
and the quality of the approximation is as better as the jump concerns a higher order derivative. 
We also refer to papers JU] , [S] , for the use of boundary layers in transmission strategy. 

Remark 1.2. It could be interesting to know if it is possible to obtain such a result in the 
non static case for which the Landau-Lifschitz equation is coupled with the Maxwell system of 
electromagnetic. For such a model an analysis of the boundary layer induced by the Neumann 
boundary condition on T is performed in p]. 

Remark 1.3. With the same method than the one used in this paper, it is possible to get global 
in time 0{e^) approximation for all s < ^. More precisely for any s < ^ there exists a family 
of solutions li^ G C°°([0, T^] x Q), of the equation Hl.l|) on [0,Ti;] x fl, of the equation (|1.3() on 
[0,T£] X r, e g]0, 1], with lim£^Q+ = oo, such that for all T > 0, there exists C > and eo 
such that for all e G]0,eo], there holds — u^\\l^([o^t]xQ) ^ Ce^ . 

2 Asymptotic expansion 

Let us fix some notations. We will use the letter S to denote the Schwarz space of rapidly 
decreasing functions. We define the boundary layer profile spaces 

Af±{T) :=F°°([0,r] xn,S{M.^)). 

Since we will need an equation of the boundary T, we fix once for all a function <I> G C°°(M^,M) 
and we assume that r2 = {$>0},r = {<I' = 0} and |V^'(x)| = 1 in an open neighborhood 
Vr of r 2. Let us also fix a function ^ G C°°(R3,M) such that S = = 0} and such that 
|V^'(x)| = 1 in an open neighborhood Ve of S ^. We assume that the neighborhoods Vr and 

^Hence for a; G O n Vr; ^{x) = dist{x,r). 
^Hence for a; G O n Ve: i>ix) = dist{x,T.). 
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Vs have been fixed small enough in order that Vr n Vs = 0- We will denote := n {'J' > 0} 
and fl- := 57 n < 0}. We consider a C°° unit vector field which coincides on Vr with 
-Vx^ ■ Va: and on Vs with -Vx^ ■ V^;. 



'J'<0 




*<0 



r 'i> = o 



In the easier case where is continuous across the hypersurface S, paper jlj shows the 
existence of solutions u^, e g]0, 1], of the equation in Q, of the equation on F, of the 
form 

where the function it is in M+{oo) and satisfies il{t,x,z) = for x ^ Vr- The function il 
describes a boundary layer which appears near the boundary to compensate the lost of the 
Neumann condition from the complete model (|l.l|) - (|1.3j) - (|1.4|) to the limit model (|1.6j) (e = 0). 
The amplitude of this boundary is weak and its behaviour is linear. For sake of completeness 
we will state this in section ?I7I\ The functions can be seen as remainders. 

Here since we deal with a ground state which is discontinuous across the hypersurface S, 
we look for solutions u^, e e]0, 1], of the equation in Q, of the equation ()1.3|) on T, of the 
form 

u'{t,x) :=U{t,x,^^^) + e(il{t,x,^^) + w^t,x)y (2.1) 

The function U describes a large amplitude internal layer profile i.e. a sharp transition in the 
neighborhood of the hypersurface S of width ~ e. More precisely the function U is C°° and 
satisfies 

lim U{t,x,y) = u^{t,x) for x G Vs n n± (2.2) 

J/— >±oo 

U{t,x,y) = u°{t,x) for X ^ Vs and y E M (2.3) 

The profile il, as we have already said it above, was constructed in |^. The functions can still 
be seen as remainders. Let us explain this time more precisely what we mean by remainders. 
Let us fix a finite set of smooth vectors fields Tq = {Zi{x\ dx)] i = 1, • " " tangent to 

the surfaces F and S (that is satisfying Zi{x]dx)^ = on F and Zi{x;dx)'^ = on S, for all 
i G {!,••• , A*}), and generating the algebra of smooth vector fields tangent to F U S. These 
vector fields can be viewed as vector fields on M'^ tangent to M x F and to M x S. By adding 
the vector field dt to the family, one gets the set T := {dt} U Tq which generates the set of 
smooth vector fields in tangent to (R x F) U (M x S). We denote Zq := dt- For all multi-index 
a G we note Z"' = d^°Z°^ . ■ ■ ■ .Z^'^ , with a = (ao, ai, • • • , a/i)- Let us introduce the usual 
norm: 

II^IU '■= X] \2'°'u\\l'2Q0^T[xQ), 

|o|<m,agNi+'' 
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and note i?™(]0, r[xf2) the space of u G L'^QO, T[x^) such that this norm is finite. We introduce 
the set E of the family (w'')o<£^i) of functions in L^(\0,T[xQ) such that for all m S N, there 
exists Sq > such that 

SUPo<e^eg (I I W ||m 

+ Ik^n W^IU + e(||w^||oo + II^W^Iloo + Ik^nW^IU)) < OO. (2.4) 

In fact Theorem 11.21 is the straightforward consequence of the following result. 

Theorem 2.1. Let u° G C°°{R,H^{Q)) be a solution of ^T^. There exist T > 0, a profile U 
in C°°((0, T) X $7 X M) which satisfies ^2.S^ — l\2.cl^ and a family (w^) in E such that the function 
given by the formula \2.1]) are solutions in C°° of the equation (| j. j|) on [0, T] x Q, of the 
equation on [0,T] x T. 

Theorem 12.11 exhibits some large variations solutions of the Landau-Lifschitz equations as 
the exchange coefficient tends to zero, by means of the asymptotic expansions (|2.1() . The 
remainder of the paper is devoted to the proof of Theorem 12.11 As in (4|, since the magnetic 
moment u is unimodular, the equation ()1.1|) is equivalent for smooth solutions to the following 
one: 

/:^(n^9)n^ = F(^i^e5^n^7^(n^)) (2.5) 

where we have noted 

C%v,d) ■.= dt-e^A^-e\ AA^, 

and 

F(n, V, H) := \V\^u + u A H - u A {u A H), 

for allu £R^, V £ M{Mp,R^), H £R^. From now on we will deal with equation (|2.5|) rather 
than Hl.l|) . We will proceed in three steps. In subsection 12.11 we will define the profile U as a 
local in time solution of a pair of nonlinear equations in Q x R± coupled by some transmissions 
conditions on {y = 0}. In subsection 12.21 we will recall the results of 1^ about the profile it. In 
subsection 12 . 31 we will prove the existence of some remainders till the lifetime T of the profile 
U. Eventually we will show that the remainders satisfy the uniform estimates uniform (|2.4() . 



2.1 Construction of the internal layers 

Despite that > when x S Q±- we will define U for all (x, z) E x M± since this will not 

cause any additional difficulty. An Uryshon argument yields the existence of two functions 
in iJ°°((0, cxd) x O) such that uf^. = for all x G il± U (il^ — Vs)- We look for a viscous internal 
layer profile lA of the form 

r u\{t,x)+U+{t,x,y) ify>0, 
U{t,x,y):={ (2.6) 
y vP_{t,x) +U-{t,x,y) ify<0. 

The functions U± are in N±{T). These functions describe some internal large amplitude bound- 
ary layers, on each side of the hypersurface S. To insure that the function U is in C^((0, T) x 
n X M) it is necessary to impose the transmission conditions: 



U+-U- = -u% + uO 
dyU+ - dyU- = 



when (t, X, y) G (0, T) x x {0}. (2.7) 
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In Theorem 12.21 we will define the profiles hl± as local solutions of nonlinear equations in 
X R-t coupled by some transmissions conditions on {y = 0}. Let us look for convenient 
equations. We will plug the functions u^''' defined by u'^'^{t,x) := U{t,x, instead of in 
(|2.5j) . In general it is not possible to verify (|2.5j) but we will try to choose the functions l{± such 
that the error term is as small as possible. Let us begin to look at the left side of (|2.5j) . With 
((221) we get in L°° 

C^u^'^, d)u''^ = dtu\ + (l{U, dt,d^)U±^ \ + 0{e) for x G 0±, (2.8) 
where the vertical bar | means that y is evaluated in y = and 

L{u,dt,d^) :=dt-d^-UAdl 

We now turn to the right side of 1)2. 5(1 . We first look at the action of TC on the family u^'^: 

niu''^) = niul) - iU±.n)\ n + 0{e). 

Then 

Y[u\ed^u',n{u^)) ■.= Y{u\,Q,n{u^±)) + F±{U±,dyU±)\+0{e) for x G !^±, (2.9) 

with for all [/ G M^, y G 7W(M^M3), 

F±{U,V) := IVI"^ {u^ + U) +U AHiul) - {U.n){ul + U) An 

+U A (n^ + U)A (W(u^) - {U.n)n) + A {U A (W(u^) - {U.n)n)) 
-{U.n)ul A (u^ An) 

Thanks to (|2.8j) and (|2.9j) we get by looking at the terms at order 

dtul + L{U,dt,d^y)l(± = F{ul,0,n{ul)) + F±(l(±,dyU±). 

Since for x G the functions satisfies ()1.6() we could simplify and we get the nonlinear 
equations 

L(4 + U±,dt, dl)U± = F±(U±,dyU±). (2.10) 

The equations 1)2. 1U() are parabolic with respect to t, y, the variable x can be seen as a parameter. 
The following theorem claims that it is possible to find some solutions U± G A/±(T) of these 
equations even for all a; G ^2. 

Theorem 2.2. There exists T > and there exist some functions U± G A/±(T) which verify the 
equations (^2.10^ when {t,x,y) G (0,T) xilxM-|- and the transmission conditions ^2^. Moreover 
precisely for all x ^ Vs and y G M± there holds U±{t, x,y) = 0. 

Proof. We will proceed in four steps. 

Step 1. We begin to reduce the problem to homogeneous boundary conditions. 
We introduce the functions V± and h(± given by the formula 

V±{t,x,y) := (l-_)4(t,x) + — nO(t,x), 
W±{t,x,y) := U±{t,x,y)±^{ul{t,x)-u^_{t,x))e'fy. 
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Thus the transmission conditions (|2.7|) reads: 

9,w+ -^wl = } e (o,r) X X {o}. (2.11) 

Moreover the equations ^UU^ - ^n} read for {t,x,y) e {0,T) x n x M±: 

L{V± + W±,dt,d^y)W± = F±{t,x,y,W±,dyW±), (2.12) 

where F± are C°° functions such that the functions F±{t,x,y, 0,0) are rapidly decreasing with 
respect to y. 

Step 2. We prove the existence of compatible initial data. 

Let us to explain why the initial values Wo,+ must satisfy some compatibility conditions at 
the corner {t = y = 0} are required in order to obtain smooth solutions W-t of the problem 
(|TT^ - (fTTT|) with W±|i=o := Wo,±. We start with the condition of order 0. Set t = in the 
transmission conditions H2.11() to see that Wo,+ must satisfy the relation 

Now, for each k ^ 1, apply the derivative d!^ to the transmission conditions (|2.7j) . We get 

8,efw; : ia^wl = } - e (0. ^) X " X {»}■ 

Now remark that, by iteration, we can extirpate df'W± by the interior equations ()2.1U|) in terms 
of derivatives with respect to y. More precisely there exists some smooth functions such that 
di'W± = C^{{dyW ±)i^2k) ■ Thus the following kth order compatibility condition must hold: 

arZn4^iLT \ \ Ak ^ht \ \ n f when (x, y) G S2 X |0}. (2.14) 
dyCli{dl^W±)i^2k)-dyC^iidl^W±)i^2k) = 0. j 

Lemma 2.1. There exist some initial values Wo,± in H^{^},S{M±)) which satisfy the relation 
IF71) and fO^ for all k ^ 1. 



Proof. As we will follows the method of we only sketch the proof for sake of completeness. 
We start by analyzing more accurately the compatibility conditions and more especially the way 
the functions depend on the derivatives with respect to y. Indeed they exists some functions 
such that 

ci((4w±),^2fc) = c'^((5[w±),^2fc-i) + {Of + (y± + w±) A af )w±. 

Since given two functions w|^'' in H°°{Q) the applications 

W± ^ W± + {V± + Wi°^) A W± 

are two automorphisms of H^{Q) and an iteration, we deduce by iteration that there exists a 
family {W^^^)km in i?°°(J^) such that 

Cl{{dyWl\^2k) - CliidyWi^^^k) = 0, 
dyC'l{{dyW'i^)l^2k) - dyC'l{{dyW'i^)l^2k) = 0. 

We end the proof by a classical Borel argument. □ 
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As a consequence, we will assume in the rest of the proof that the functions W.^ vanish for 
t ^ 0. 

Step 3. We look for linear estimates. 

In order to use an iterative scheme, we look at the linear problem 

L{W±,dt, d'^)W± = /± when (t, x, y) e (0, T) x S7 x M±, (2.15) 



W+ - W_ = 0, 

dyW^ - dyW^ = 



when G (0,r) X X {0}. (2.16) 



For all real A ^ 1, the space L^((0,r) x^x M.±) is endowed with the scalar product associated 
to the Euclidean norm 

l|W±||o,A,T := ||e~^* W^±||L2((0,T)xnxR±) 

In order to avoid heavy notations, we will denote W := (W_|_,W_), / := (/+,/_) and SZU := 
(im+,2n_). We endow the space L2((0,r) x x M+) x L^((0,r) x J7 x M_) with the scalar 
product associated to the Euclidean norm 

ll^llo,A,T := ||W+||+^o,A,T + ||W_||__o,A,T- 
For m G N, we introduce the following weighted norms: 

\\W\U,x,t:= Y1 P"x^IIo,a,t, and |iyU,A,T := ^1 lo,A,T, 

where a := (ao, as) G and := d^°d^'d^'d^\ 
Proposition 2.1. Let R> 0. IfW± verify the following estimates 

l|2n+||ijp((o,T)xQxR+) + l|2n+llLjp((0,T)xnxK+) + |2B|m,A,T < R-, 
and the following boundary conditions 



2IJ+ - 2n_ = 0, 
dyW+ - dyW- = 



when (i,x,y) G (0,r) X 1^ X {0}, (2.17) 



then there exist \m > and for all A; G N, fj,k,m > 0, such that for all A ^ Xm, 

\W\mXT ^ ^\f\m,X,T (2.18) 

and for all fi ^ fJ-k,mj 

k 

|y' W\m,X,T ^ ^ ^ |yi (2.19) 

Proof. We multiply the equation ()2.15p by W-t and integrate for {x,y) € Q x Mj-. Hence 

(1/2)9^/ |W±|2-Ji,±-J2,±= / /±.W± (2.20) 

JnxR± JnxR± 

where Ji,± := / W±.5^W± and J2,± := / W±.(2B± A 5^)W±. 
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Integrating by parts, we get 

Ji,± = - [ \dyW±\^ - /i,±, and J2,± = - [ W±.{dyW± A dy)W± - h,±, 
where Ii,± := f (W±.9j,W±)|y=o, and Is, ± := / ^ dy^W ±))\y=o. 



Using the boundary conditions (|2.16p and (|2.17|) . we get — = l2,+ — I2,- = 0. Taking 
that into account we add the two estimates in H2.2U() . Then we multiply by e"^"^* and integrate 
in time. By a Gronwall lemma we get that there exists c > such that for all A ^ c, 

\dyW\l^^T + m\lx,T ^c\<f,W >A,T I- (2.21) 

We go on with estimates tangential to {y = 0}. To do this we apply the derivative d"^ to 
the equations H2.15() - H2.16|) . So we get that S"^, W± verify 

L{W±,dt, dl)dl^ W± = /± when (t, x, y) e{0,T)xQx M±, (2.22) 



when {t, X, y) G (0, T) x x {0}, (2.23) 



where 



/±:=5,%/±+ E (2.24) 

|oi| + |a2| = |aMa2|<|a| 

We apply the tangential derivative df^ to the boundary conditions (|2.17|) and get 



9,",2IJ+-9,°,2IJ_ = 0, 

dydf^,mi+-dydf^,W.=^ 



when (t,rc,y) G (0,T) X J7 X {0}, (2.25) 



Using the estimate (|2.21() . we get, for all A ^ c, 

W\lx,T + W^lo,A,T ^ c| < /, 5,^, >a,t I- 

Thanks to we get 

</,9,%W^>A,T=<9,-,/,5,^,W^>A,T+ E (2.26) 

|ail + |a2l = |«l,la2|<la| 

where 1^1,^2 ■= I+,ai,a2 + ^-,01,02 ^ith 

/±,ai,a2 :=< 2IJ± A ^2 a,"^ W±, a,"^ W± >A,T . 
Using Cauchy-Schwarz inequality, we get 

I < dr^,f,d^^,W>x,T I ^ |/|0,A,T-|W^|0,A,T. 

We are going to estimate, for all 01,02 such that |ai| + |a2| = |a|, |a2| < \o(\, the term Iai,a2- 
Integrating by parts, we get /±,«i,a2 := ELi -^±,ai,Q:2' '^i*^ 

4,ai,a2 -<5,^i5,2n±A9,a,"^W±,9,",W±>A,T, 



7^2 

±,ai ,a2 
-'±,01,02 



- < d^^l 2B± A dy d^^l W±, 5,-, 5,W± >A,T, 

T « {(5i"i2n± A9,5,"^W±)}|,=o,K,9,W±}|,=o »a,t, 
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where << .,. >>a,t denotes the scalar product of L^((0, T) x Q) associated to the mesure 
e~^^dtdx. Thanks to the boundary conditions we get I+^ai,a2 ~ -^-,ai,a2 = ^• 

Using Cauchy-Schwarz inequahty, we get 

\lia,,a,\ ^ \dr,'^dyW±Adyd^^lW±\o,X,T.m±\U,X,T, 
\lla„a,\ ^ |a,"i2IJ±A5,a,^^W±|o,A,T.P,W±|U,A,T, 

Using Garghardo-Nirenberg inequalities, we get 

l-^i.ai.aal ^ c{\\dyW±\\m,\,T-\\^±\\Lip + \\^±\\Lip-\\dy^±\\m,\,T)-\\^±\\m,X,T, 

m—l,X,T- m—l,X,T m,\,T- 

Hence we get 

We deduce that there exists Am > such that for all A ^ Am, there holds ||Vl^||mAT ^ 

^\\f\\m,X,T. 

To prove the estimates ()2.18p . it remains to get normal estimates. The cases 04 = or 1 
are already treated in the tangential estimates. If 04 ^ 2, we proceed by iteration, extirpating 
dy'W± from the equations. 

It remains to get the estimates ()2.19|) . First we notice that for p ^ 1 the function W-t 
verify the initial boundary value problem 

L{W±,dt, 9^)W^^ = /l'^ when (t, x, y) € (0, T) x x M±, 



w!f] - wi?] = 



■ when (t, X, y) G (0, T) x Jl x {0}, 
= when (t, x, y) G {0} x 17 x M±, 



where 



p-i 
i=o 

where the and the g| are in N. Thus we prove, by iteration on p and thanks to the inequality 
H2.18() . the estimate 

p 

j=0 

which implies the estimate ()2.19() . 

□ 

Step 4. We use an iterative scheme. 



We define the iterative scheme (W!^)i/gN by setting Wj. equal to zero and, by iteration, 
len is defined, we take W!^"^^ as solution of 

L{V± + W^,dt,dl)W^^ = F{t,x,y,Wl,dyW'^) when {t,x,y) G (0,oo) x 17 x M±, 



Wl+^ - Wi^^ = 0, 



dyWl+^ - 5yW^+i = 



when it,x,y) G (0,T) x 17 x {0}, 
W;^;*"^ = when {t, x, y) G {0} x 17 x ]R±. 
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Thanks to the Unear estimates, to a Sobolev embedding and to some Garghardo-Nirenberg 
inequahties, we show that the iterative scheme (W5j.)jygi^ converge, when v — > +00 toward some 
solutions W-i- G A/±(T) of the problem ()2.12|) - ()2.11l) . By going back to the original problem 
H2.10() - H2.7|) . the first sentence of Theorem l2.2l is now proved. When x ^ Ve, the function — ti*^ 
in the right hand side of (|2.10|) vanishes and so do the functions IA±. 

□ 

Remark 2.1. Notice that the possibility of a blow-up can be controlled with Lipschitz norm 
in a very classical way. However we do not know whether the solutions lA actually blow-up or 
globally exist. 

2.2 Construction of IX 

In this section we define the boundary layer profile 11 as a solution of a linear boundary value 
problem. Let us recall that this function describes a boundary layer which appears near the 
boundary to compensate the lost of the Neumann condition from the complete model Hl.l() - 
(|1.3|) - (|1.4|) to the limit model (|1.6j) (e = 0). Such a boundary layer was already mentioned in 
paper RP. Let G be a C°° function on such that = 1 in a neighborhood Wr of F such that 
Wr CC Vr and B = in 17 - Wr. 

Theorem 2.3. There exists il G A/'+(T) which verifies 

L{u^,dt,dl)^ = -(ll.n)n° An + llAW(u°) 

+H A (n° A n{u^)) - (U.n)n° A (n° A n) + n° A (H A W(n°), 

when {t,x,z) e (0,r) x f7 x M+, 

d^il = e{x)dnu° when {t, x, z) e {0,T) x Q x {0}. (2.27) 

Moreover there holds x,z) = for x ^Vy,- 

Proof. Proceeding as in the proof of Theorem 12.21 we prove the existence of compatible initial 
data. Then we follow the proof of Proposition 4.2 of □ 

2.3 Construction of 

In this section, we look at the remainder w^. We will proceed in four steps. First in section 
12.3.11 we will reduce the initial problem Hl.l|) - (|1.3() - H1.4|) for the unknown to a problem for 
w^. Indeed in order to get a homogeneous boundary problem, we will add a corrector to 
w*" and rather work with the resulting term . Some Borel classical arguments will insure 
the existence of convenient initial data for the resulting reduced problem which means that 
compatibility conditions either on T and on S are satisfied. We will prove that the solutions of 
this nonlinear problems exist not only for a common non trivial time, in fact even till the lifetime 
T of the profiles li. Moreover these solutions satisfy some estimates uniform with respect to e. 
The method lies on a simple Picard iterative scheme (cf. section r2. 3. 2|) and on linear estimates 
(cf. section [2.3.3)1 . More precisely we will use L'^-type conormal estimates of only the two first 
normal derivatives, and some Lipschitz estimates. A few carefulness reveals that the presence 
of the operator 7i does not cause any loss of factor e or any loss of derivatives. 
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2.3.1 A reduced problem 

Since we look for solutions of (|1.1|) - (|1..S)) - (|1.4|) of the form (|2.H) where the functions 

have been constructed above, we look for a problem in term of the remainder w^. In fact, in 
order to get a homogeneous boundary problem, we choose a function p[t, x) G H°° such that 

5n/o|r = -anH(t,2;,0)|r. (2.28) 

and will look for remainders of the form w'' = p + w^ . Let us explain why. On the boundary 
r, the function satisfies: 

9na^|r = ea„il(t,x,0)|r, (2.29) 

Hence in general does not satisfy the homogeneous Neumann boundary condition on T. We 
define the function := + ep. Thus we look for solutions of Hl.l() - H1.3() - H1.4() of the form 
tt"^ = + e = a!^ -\- e . Combine ^\.\^ , (|2.28|) and (|2.29|) to find a homogeneous Neumann 
boundary condition on F for w^: 

d^w' = on]o,r[xr. (2.30) 

We now look for an equation on the unknown w^. The function belongs to C^((0, T) x 0) 
and to H^{^}). Moreover, d^ satisfies the equation 

£(d^5)d^ = F(d^ea^d^?^(d^)) +er^ (2.31) 

where the family {r^)e lies in the set E (defined above Theorem 12.1(1 . The system for the 
unknown w'^{t,x) writes 

C{¥ + ew^,d)w^ = K{e, d^e9^d^ W(d^), ed^w^,n{w')) + r= in ]0, T[xn (2.32) 

where ii' is a smooth function of its arguments. Let us use more concise notations, and note 

:= (d^ea^d^?^(d")) and W := {w' , ed^w' ,n{w')) . (2.33) 

Then, the Taylor formula shows that the function K has the following form: 

K{e, A^, = G{e, A^,eW^)W^ 

where G depends smoothly on its arguments (including e), which will be useful in the sequel. 

Following 13_ there exist a family {wf^:^^)e of compatible initial conditions for the problem 
(j2.32j) - (|2.30|) which verifies suitable uniform estimates with respect to e. We choose such a 
family. 

2.3.2 The iterative scheme 

We want to solve the problem (|2.32() . (j2.3()|) . We use a simple Picard(-Banach-Caccioppoli) 
iterative scheme defining a sequence w^'^ which will converge to the solution of the problem. 
For clarity, we adopt the following more concise notations 

A' := (d%ea^d%W(d^)) and W'^"" := {w''"' ,ed^w'''' ,n{w'''')) . 
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With these notations, the iterative scheme writes 

C{cf + ew^-", d)w''''+^ = f''' in ]0, T[xn (2.34) 

where 

f'" := G(e, A',eW'''')W'''' + (2.35) 
This equation is coupled with the initial and boundary conditions: 

dnw^'^^^ = Oon]0,T[xr (2.36) 
w^'-'+^^o = <if (2.37) 
The iterative scheme is initialized with w'''^{t,x) := wf^^^{x). 

2.3.3 Estimates for a linear parabolic system 

Consider the linear problem 

£(a^ + eb,9)u = f on]0,T[xn (2.38) 
dnU = Oon]0,r[xr, (2.39) 

We endow the space H^^{]0,T[xil.) with the usual weighted norm with A > 1: 

In order to estimate the initial data, we introduce the similar norms built with the set Tq instead 
of T, integrating on Q instead of [0, T] x Q: 

|-,| \ ^ |o| II 

|o|<m,ao=0,oGNl+f' 

We will use the following classical Gagliardo-Moser-Nirenberg estimates for conormal derivatives 
(see 0). 

Lemma 2.2. Let m G N. There is Cm > such that, for any ai, . . . ,ak € HJJ^(\0,T[xQ) n 

L°°{]0,T[xQ), for all multi-index qi G N^+^ . . . , G N^+^ with \ai\ H h \ak\ < m, for all 

A > 1.- 

||Z"iai...^"'=afc||o,A < c„, Yl {hj\\m,xll\\ai\\oo)- (2.40) 

The following proposition gives some e-conormal estimates for the two first normal derivatives 
of the solutions of the problem 

Proposition 2.2. Let R > be an arbitrary constant and m ^ 3. There exist Cm(i?) > and 
Am > such that for a fixed constant large enough, depending only on the choices of the vector 
fields Zj, the following holds true. Assume that 

^( l|b||oo+ Il^jblloo + Ik^^blloo ) < (2.41) 

then, for all A > Am, the following estimates hold: 

lk5a;U||m,A + A||u||m,A < C'm(^) [ ||/|U,A+ ^m,A(u) 

+ e ( ||e(9xb||m,A + ||b||m,A ) ( ||u||oo + |k<9a;U||oo ) ] , 
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where 

Va(u) := ^ \{dtu)\t=o\m-k,\- 

0<k<m 

and 

Wi^dn) u\\m,X < Cm{R)[ \\f\\ m+l,X + e||b|U+l,A (||u||oo + 

+e||eanu||„+i^A + e^||u|U+2,A 
Proof. Step 1. Let us note v := e""^* u, which satisfies 



(2.43) 



£(a^pp + eb,9)v + Av = e"^*/ on ]0, T[x!^ (2.44) 
(9nv = Oon]0,r[xr. (2.45) 
V = wf^^^ on t = 0. (2.46) 

Let us note ||.||i2 the L norm in [0, T] x VL, and |.|^2 the L norm in VL. Multiplying (|2^ by 
V and integrating on ]0,T[xf] gives the following estimate, integrating by parts the e^/S.x with 
Green's formula in Vl: 

^'l|V..v||i2+A||v||i2 <2 |((e-^*/,v))i2| + |v(0)|i2, (2.47) 
for all A > Ao if Aq is fixed large enough, and for all e > 0. In terms of u it writes 

e^\\VM\lx + A||u||2^ < 2 |((/,u))^2j + |u(0)|i2, (2.48) 

where L| is the Hilbert space L^(]0, T[xr2, d/i) with the measure dpL := e~^^*dtdx. 

Using now the Cauchy-Schwarz inequality in the right hand side, and absorbing in the left 
hand side the term in ||f ||^2 yields the desired estimate for m = and some constant cq > 0. 

Step 2. We show the inequality by induction on m. Assume it for m — 1. We apply a tan- 
gential operator with fields Zi to the system, and \a\ = m. The function Z^u satisfies 
the same boundary conditions. The estimate (|2.48|) gives, for A > Aq: 

e2||V,Z"u||2, + X\\Z^u\\l, <2|((e-^*Z"/+ [(a^.pp + eb)e2A„Z-] Au,Z"u))^2|. (2.49) 

where [., .] denotes the commutator. Using Cauchy-Schwarz inequality and 2ab < 2A^^ a^+A6^/2 
yields: 

^2||v7 Q-a^llZ I '^||-7a„||2 ^ ^ \\ -Xt ■^ra f\\2 

e llVxZ U\\r2 + —\\Z U\\r2 < — e Z J\\r2 

II ^ UL 2 A ^ (2.50) 

+ 2 KiMpp + eh)e'Ax, A u, Z°u))^2 1 . 

We need to control the second term in the right hand side of (|2.5U|) . The commutator [ a^^p e^A^;, Z'^ 
writes as a finite sum 

£2 ^ a^(t,x)Z'3 + e ^ 6^(t,x)e5„ZT+ ^ 4it,x){ed^)^Z^ (2.51) 

|/3|<m+l |7|<m \S\<m-l 

where the coefficients a^^ , bf^, c| are bounded functions satisfying 

sup ||e«9na^||Lcx>(n) + ||e(9n6^||L°°(f7) + Ik^ncf ||loo(q) < oo (2.52) 

£G]0,1] 
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for all P,j,5, because (|2.5'2j) holds clearly if we replace L°°{Q) by L°°{Q^) or by L°°(r2_), and 
because a^^p is in H^{Q) for all e > 0. Hence we are led to control the corresponding three sort 
of terms: 

e2(( a^Z^u, Z-u e(( b^edn)Z"'u, Z^u (( 4{edr.fz'u, Z^u (2.53) 

where \f3\ < m + 1, \j\ < m, \6\ < m — 1. The first two terms in H2.53() are simply controlled by 
(5||eVa;u||^ ^ + Cs6^^ ll^llmA ^'^^ ^ arbitrarily small, and Cs being a constant depending on 6, 
but independent of e. For the third term one uses an integration by parts (by Green's formula) 
of the field to show that this term writes as a sum of terms of the form 

e^~^~^'{{{ednyz'u, {edny'z"vL))^2 

where \6\ < m — 1, j,j' £ {0,1}, and is a bounded function (uniformly in e) since all the 
boundary terms terms vanishes: dnZ°'u\g^ = 0, for all a G M'^. It follows that the third term in 
(|2.53|) is controlled by CA~^||eVj;u||^ + C||u||^;^ for a constant C independent of e, and all 
A > 1. Hence, by choosing a 6 > arbitrarily small, and Ai > large enough, there holds 

I (( KppS^^ccZ''] An, Z"u))l2 I < 5||eV,u||^^;, + c^||u||^^;, 

for all A > Ai, and for all e g]0, 1], with a constant Cm independent of e. 
We need now to estimate the term 

(([ebe^A,,^"] Au,Z"u))i2. (2.54) 



The commutator [ h e"^ Ax, Z'^] writes as a finite sum 

E a,M2'h)Zf^' 

\l3\<m,\f3'\<m+l,\l3\ + \f3'\<m+2 

+ e Yl 6^,y(Z^b)(e9„)Z^' 

|7|<m,|7'|<m,|7| + |7'|<m+l 

+ Yl cs,s'{Z'h){ednfz'' 

\5\<m,\5'\<m-l,\S\ + \S'\<m 

where a/3,/3', ^7,7', cs^s' are smooth fonctions on 0,. Hence to control the term ()2.54|) we are led to 
estimate tri-linear terms in (b, u, u) of the following form (where dfi := e~'^^^dtdx): 

f ap^(3'2:%.Z'^'ui.Z''uj dn, |/3| < m, |/3'| < m + 1, + < m + 2 (2.55) 

J]o,T[xn 

Z^h . ednZ^'u, . Z"uj dfi, I7I < m, |7'| < m, \j\ + < m + 1 (2.56) 



T\xn 



[ cs^s' ^^b. {ednfZ^'ui .Z^Uj dfi, \6\ < m, \6'\ <m-l,\6\ + \6'\ < m, (2.57) 

where the Uj are the components of the vector u. Let us treat the term ()2.57() . By the green 
formula, the integral can be written as a sum of integrals of the form 

[ cs,s' Z^ednh.Z^'ednUi.Z'^Uj dfi (2.58) 

J]o,Tixn 

[ cs^s' 2:^'b-2^'£dnU,.Z''ednUj dfi, (2.59) 

e [ ds^s' Z^h.Z^'ed^Ui.Z'^Uj dfi, (2.60) 
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and other terms involving lower order derivatives easy to control. The term ()2.58|) is controlled 
by 

which is bounded by using the Gagliargo-Nirenberg-Moser estimate by 

c( \\£dnh\\m,\\ednU.\oo + \\edn\l\\m,\\l^dnh\oo ) \\'H-\\m,\ 

and hence by 

C(l + R){ ||e«9nb||m,A|e<9nU|oo + ||e9nU||r„,A) || 

For the term H2.59() there are two cases. The first case is when 6 = 0. In that case the integral 
is bounded by 

The second case is when \5\ > 1. In that case we write Z^h = Z^h and apply the Gagliardo- 
Nirenberg-Moser inequality with Zh in L°°. The term in bounded by 

c( \\Zh\\m-l,x\ednU\oo + \\ednU.\\m-l,\\Zh\oo ) \\£dnU\\m,X 

and hence by 

c ||b||m,A|e9nu|oo|k9nu||m_A + cRX~^ \\ednu\\l^^x 

The next terms like H'2.6U() are easier to treat in the same way, and are bounded by the same 
terms. The term (|2.57|) was the more delicate to estimate. The terms (|2.56|) and (|2.55|) are 
simpler and can be treated in a similar way. Replacing in the right hand side of 1)2. 49() and 
summing over all the possible operators -Z" gives the desired estimate, and the proposition is 
proved. □ 



2.3.4 Iteration 

Now classical arguments show the convergence of the iterative scheme if e g]0, Eq] and eq is small 
enough. We describe the main lines (see jl3j). Let us fix an integer m > 4, and note 

R:=l+ sup {£( ||z/;^'0||oo+ V H^jW'^'loo + Ik5,w;^'0||oo )}• 

Proposition 2.3. Let be given A > 1. Then there exists h > 1 such that for Eq > small 
enough, for all G N, for all E G]0,eo], there hold 



k^'lloo 



+ \\ZX 'l\oo + \\ed^w'''' Woe <Re-^ (2.61) 

0<j<M 



and 

lk''1U,A + IkSnio^'HU.A < h. (2.62) 

Proof. For h large enough, the inequalities 1)2. 61() and 1)2. 62|) are satisfied for 1^ = 0. Now suppose 
that w^'" satisfies (tTTm . TUrlh . We want to prove that ■w^-'^+i also satisfies (tTTm . (tTO . 
The proposition 12.21 gives a constant Cm{R) and the inequality (|2.42|) holds with u = w^''^'^^, 
b = w'^''^, and / = /^''^ defined in (|2.35|) . In order to control the right hand side of l|2.H4j) . 
we need a control of ||'H(i(;^''^)||oo and of ||'H(if^''')||m,A) which is a consequence of the following 
lemma. 
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Lemma 2.3. Let m € N. There exists c > such that for all A > 1, 

\\n{v)\\m,x + |k5nW(t;)|U_i,A < c{\\v\U^x + \\edr,v\\m^i^x) ■ (2.63) 

Proof. We note E{d) := (div ,curl ) the operator from [5'(M^)]^ to We denote by 

E~^{d) the inverse operator. Then u = E~^{d)f , is defined by u{£) = {a{£)i — i ^ KO) 

where /(^ = {a{0,HO) G M x M^. Thus n(0 = M{^)f{C), where M(^) is a 3 x 4 matrix whose 
entries are rational functions of homogeneous of degree —1. Let us fix x G C^(M'^,M) such 
that x(0 = when |^| < 1 and x{C) = 1 when |^| > 2, and call P{D) and R{D) the operators 
from [5'(M3)]4 [5'(r3)]3 defined by P{D)f := J^~^{xMf) and R{D)f := (1 - x)f) 

where J-~^ means the inverse Fourier transform. In the sequel we will simply note 5'(M'^) 
and L'^{n) instead of [^'(M^)]'^ and [L'^i^)]' , meaning that we talk about the components of 
the vector valued functions, the (finite) number of components being understood. We have 
E~^{d) = P{D) + R{D). The operator P{D) is a special case of classical pseudo-differential 
operator of class S^-^ o(^'^ ^ elliptic, and R{D) is an infinitely smoothing operator of class 

5{;^(m3xm3). 

Let us now take into account the t coordinate. Let us note Q =]0, r[xr2, F =]0, r[xF and 
S =]0,T[xE. We extend the actions of P and R to the spaces of functions or distributions 
which depend also on t like L'^{d) or C([0,T],5'(M^)), by considering t as a parameter so that 
Pu{t,x) := P{D)u{t, .){x). Let v G H^iQiR'^) such that dnV G H^-\n;R'^). Then n{v) = u^^ 
where u G L^i[0,T] x M^) is defined by E{d)u = E{d)v + (f|p.n) 6f,, where the notation V 
means the extension of F by to [0, T] x (or to R^, depending on the context). 

Let us note / := E{d)v, which is in H^~^{]0,T[xM.^), and g = (f.n)|r. This trace is 
well defined since by assumption v G H^{Cl), and using local coordinates patches one sees 
that g G i?™'~2^r), the usual Sobolev spaces. The operator P{D) satisfies the transmission 
property (introduced by Boutet de Monvel jj, [2]) on Q and on \ O because its symbol is a 
rational function of ^, which is a sufficient condition to satisfy the transmission condition. The 
transmission property has been also studied and used by Grubb, and we also refer to papers [HI 
and j2|. To avoid many repetitions, we will note in what follows Qi := Q and ^2 = \ Since 
P{D) is elliptic of order 1, the transmission property implies (see and [7]) that if t> G H^{Q) 
then for j = 1,2, (^P{D)v) G H'^~^^{^j). Let us note note uq) = u^^,, for j = 1,2, so that 

■H{v) = U(i) G L^{n). Using the notations of P, [S], [7], 

= {E-\D)v)^^^^^ = P{D)(^^f + Ki.'\g) + R{D)^^^ E{d)v, (2.64) 

where P{D)^^^f = {P{D)f)^^^^^, R{D)^i^v = {R{D)v)^^^^^ and where K^^\g) = {P{D){g ® 
) "Poisson operator": 

K^^^ : H'{T) ^ i?^+i/2(rj(j')), (2.65) 

(linear continuous), extended to functions depending on t as a parameter. (See theorems 2.4 
and 2.5 of 13). 

Let us now prove the lemma. First of all, d^Tl(v) = Tl{dY^v) is in L^{Cl) because (9™u G 
L^(i7) and H acts on L^{tl). It is also easy to show that d^~^H{v) G H^{0,): by assumption, 
for any t G [0,r], -^t;(t, .) G H^Q), hence n{d'^'~^v){t, .) G H\n) because ^^-^^(t, .) is 
piecewise-i/^ and because of the properties of H. Hence dxd^~^7i{v) G L?'{^) and since we 
aheady know that dY'Hiv) G L'^{^) we have proved that d'^-^n{v) G H^{Vt). 
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Let us show now that Zjd'^~'^n{v) G H^{n) for j = 1, . . . We have E{d)u = f + g® 
5y- Since E{d) is ehiptic (as an operator in 5'(]R^), but not in we can express the 

normal derivatives of u in term of tangential derivatives and of E{d)u, and this implies that the 
commutator [E{d), Zj]u writes 



[E{d),Zj]u = 2^ AjZjU + Aof + Bg(S)5r (2.66) 
1 

where Aj , B are matrices with entries (depending on the fields Zj ) . It follows that 

E{d)ZjU = M^Z'^f + N^iZ'^g) ® 6r 

\a\<l \a\<l 



with C^(M'^) matrices Ma,Na, and applying gives: 

E{d)Zjdr~^u = Y Ma-Z^^d^-^f + Y N^-iZ'^dr'^g) 5r (2.67) 

|a|<l |o|<l 



Now Z'^dJ^ ^/ G L^(il), because / = E{d)v, and the transmission property implies that 
for every t G [0,r], the function P(L»)(j) ( "V ) (*, ■) is in iJi(OQ)). This implies that 
d^P{D)^^^ ( - V ) g ^2(j^) gij^pg already know that dtP{D)^^^ ( " V ) G L'^{^) 
from the previous case, we deduce that for j = 1,2 the functions P{D)^^^ (^Z^d]^^^ f ) is 
in i/^(i7(j)). Concerning the boundary term in (|2.67j) . since g G H^~^{T) we know that 
Z°'d1^~'^g G and the property implies that, for ah t G [0,r], the functions 

Ki.^\z''dr~^g){t,.) is in ff^(r2(j)). By the same way as before we deduce that for j = 1,2 
the functions K^^^Z'^d^'^g) is in H^Q^j)). Now, applying E{d)-^ = P{D) + R{D) to the 
equation (|2.67p gives Zjd^~'^U(^j) G ff^(f2(j)) as claimed. Then, the proof can be continued by 
induction in the same way. □ 

The lemma 1231 together with the Gagliardo-Nirenberg-Moser estimates and the induction 
assumption, implies that (like the majoration of the term (5.25) in paper !■{ ) : 

||r'1U,A < ciR)i\\w''''\\m,x + lk9,^''1U,A) < c(i?)p(A). (2.68) 
Hence, the proposition 12 . 21 implies that 

\\ed^w''''+'\\m,X + X\\w''''-^'\\m,X < Crr,{R)[X'' c{R) p{X) 

We now use the following Sobolev inequalities: 

e^/^||u||oo < e"\||u|U,A + ||e5nu|U,A), 

£'/'lk5„u|U < e''\\\ednU\\m,X+\\{ednfu\\m,x). 

By taking A large enough et e > small enough the inequality 1)2. 61(1 is also satisfied for w''''^'^^ 
and the proof by induction is complete. □ 

Now by extracting a convergent subsequence it is a classical argument to show the conver- 
gence m L2(]0,r[xl7) of w^'" to a solution of the non linear problem which satisfies the same 
estimates (|2.61j) , (|2.62|) . This concludes the proof of Theorem 12.11 □ 



(2.69) 
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